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Abstract 

From a general ansatz for the effective potential of cosmological perturbations we 
find new, exact solutions in single-scalar-field inflation: a three parameter family 
of exact inflationary solutions that encompasses all exact solutions that have 
been known previously (power-law inflation, Easther's model, and a generalised 
version of Starobinsky's solution). The main feature of this new family is that 
the spectral indices are scale dependent. 
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1 Introduction 



The observation of the first acoustic peak in the spectrum of anisotropies in the 
CMBR (Cosmic Microwave Background Radiation) ^[ leaves inflation as the 
only mechanism that can provide seeds for the observed large-scale structure of 
the Universe. The theory of inflation predicts primordial density perturbations |§ 
and primordial gravitational waves 0] from quantum fluctuations of the vacuum 
during the inflationary epoch. So far, inflationary predictions have been based 
either on the slow-roll approximation || or on power-law inflation |J. However, 
it is clear that inflation might occur in a much broader context and we should 
discriminate between generic predictions of inflation and predictions of a specific 
scenario. 

For any model of inflation, the scalar power spectrum can be expressed as a 
Taylor series around some pivot scale k* |7], || 



ln(,eP)=ln(eP*) + K-l)ln(^ ' <j " 



J 2 d In k 



A similar expression can be written for the power spectrum of tensor perturba- 
tions. The quadratic term represents the first deviation of the spectrum from a 
power-law shape and is often referred as the "running" of the spectral index n. 
No presently known exact solutions || gives rise to a "running" spectral index. 
In this letter, we present a class of exact solutions which do exhibit /^-dependent 
spectral indices. 

It is difficult to predict the coefficients of the series (|IJ) in general. One has to 
rely either on approximations or on numerical mode-by-mode integration [II], [TT], 
[T2| . The best currently available analytical method is the slow-roll approximation 
H, [L3[ where "running" only shows up at second order in the slow-roll parameters 

|8|. The precision of any approximation scheme should be controlled, this 
means that the errors should quantified. This can only be done by comparing the 
approximate power spectrum to an exact one. For inflationary models without 
"running" of the spectral indices, we have already analysed this question for the 
slow-roll approximation up to the first order by comparing to the exact solutions 



of power-law inflation |TJ], [15] . Therefore, to go beyond the first order and to test 
the accuracy of methods that estimate the "running" of the spectral index, exact 
solutions with ^-dependent spectral indices are needed. So far, only numerical 
solutions are available to perform an error analysis model by model for the second 
order of the slow-roll approximation and for models where the slow-roll conditions 



are not satisfied [[U], [TT], [Tj| . The exact solutions presented here constitute a first 
step towards a more complete analysis. 

The generic features of inflation might be studied from two sides: from the 
particle physics point of view the starting point is the inflaton potential V(ip), 
together with initial conditions for uo and up. From the cosmological side, one 
might start from the effective potential (or effective action) for the evolution of 
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cosmological perturbations, together with initial conditions for the Hubble rate 
H and its time derivative. Here we choose the latter approach, which was also 
taken in Refs. || lt\. We find a family of exact inflationary solutions for the 



mode equations of cosmological perturbations that includes all previously known 
solutions as limiting cases and that allows us to gain new insight into the generic 
properties of inflation. 

Let us now introduce some basic tools for the study of cosmological pertur- 
bations. During inflation, density (scalar) perturbations and gravitational waves 
(tensor perturbations) can be characterised by the quantities Hs,t(v) EH) which 
obey the equations (we use the notation of |L8L \U 



ii , 

PS,T ' 



k - 



y" 



A*s,t = 0, 



(2) 



where a prime denotes differentiation with respect to conformal time rj and zs,t 
are functions of the scale factor a{rf) and of its derivatives only: 



(3) 



with 7(77) = 1 — (a 2 /a' 2 )(a' /a)' , which is nothing but the slow-roll parameter 
e = —H/H 2 . Inflation occurs when 7 < 1. Eq. (0) may be viewed as a "time- 
independent" Schrodinger equation with an effective potential given by 



Us,t{v) = 



(4) 



-jfc(r/-97i) 



2k 



(5) 



The assumption that the quantum fields are initially placed in the vacuum state, 
when the mode k is subhorizon, fixes the initial conditions for the quantities /xs,t- 
They read: 

lim I^s,t(v) = =F4\/7r/pi— 

where 771 is an arbitrary initial time at the beginning of inflation. Then, an in- 
tegration of Eq. (g) allows the determination of the power spectra P{k). For 
density perturbations we choose to work with the power spectrum of the hyper- 
surface independent quantity (, which corresponds to the perturbation of the 
intrinsic curvature of a spatial, uniform density hypersurface. For gravitational 
waves we determine the power spectrum of the amplitude h. These power spectra 
are calculated according to: 



k 3 P, 



k 3 




8vr 2 





k 3 P h 



2tf 

IT 2 



(6) 



Measurements of the CMBR anisotropics and of the large-scale structure probe 
scales which were well beyond the Hubble radius at the end of inflation. Thus, 
we are interested in the modes which satisfy kj (aH) <C 1 at the end of inflation. 
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It is easy to see from Eq. (0) that the spectra are time independent in this regime 
since /ig,T oc ^s,t once the subdominant mode can be neglected. Let us also note 
that the spectrum of the quantity ( is related to the spectrum of the Bardeen 
potential today by k 3 P^ = (25/9)A; 3 P$. In agreement with Eq. ([!]) the spectral 
indices are defined by: ng — 1 = dln(/c 3 P^)/dln k and nx = d \n(k 3 Ph)/d In k. 

The behaviour of cosmological perturbations during inflation is completely 
fixed once the functions z StT (rj) are specified. The effective potential Us t T{v) 
follows from Eq. ([|) and Eq. (|2|) can be solved. Integration of equation (|3|) 
provides the scale factor. Instead of assuming specific functions ^s,t(^), one can 
also start from the effective potential ?7s,t(^) itself, such that Eq. (0) may be 
solved analytically. In this case the functions zs^iv) are known explicitly since 
Zs,t(v) = /is/rO = 0,77). 

Let us now discuss the shape of £/s,t(??)- One expects that the effective poten- 
tial possesses the following properties: lim^|^ +00 Us,t(v) ^ ^ 2 > f° r an Y wavenum- 
ber of interest, in order for the modes to be inside the horizon at the initial time rji 
(so that we can fix the normalisation from quantum- mechanical considerations); 
it seems also reasonable to assume that lim^^o ^s,t(*7) 3> k 2 which guarantees 
that the quantum fluctuations are amplified (i.e. are frozen outside the horizon). 
These requirements do not single out a unique effective potential. However, a 
general simple ansatz satisfying these conditions is given by: 

UsAv) = E (?) 



m=l 



where the c m 's should be determined for each specific model of inflation. It is 
necessary to have c m j n and cm positive. The coefficient c\ defines a characteristic 
scale kc ~ c%. The corresponding term in the series (0) dominates the other ones 
when 1 77 1 — > 00 and will therefore determine the power spectrum at very large 
scales. 

Inflationary models for which the evolution of cosmological perturbations has 
been studied so far, power-law inflation and slow-roll inflation, are chosen such 
that c m = if m 7^ 2. In this sense they are special models. In particular, they 
lead to /c-independent spectral indices. 

In this article we study the behaviour of cosmological perturbations for models 
characterised by an effective potential given in Eq. (|7|). A first step in the analysis 
of these models consists in studying the effective potential 



\r)\ \r]\ z 



The final power spectra depend on the free parameters C\ and c 2 plus one pa- 
rameter contained in z^^{r\). We demonstrate that the general solution for this 
potential is given in terms of Whittaker functions and we show explicitly how all 
the previously known cases can be recovered from this more general solution. 
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The initial conditions together with Eq. (^j) determine the region of the infla- 
ton potential explored during the evolution of the field. We always take ji < 1, 
such that an inflationary phase is guaranteed. Although all the properties of the 
cosmological perturbations can be calculated analytically, it is not possible to de- 
termine the corresponding inflaton potential analytically. This issue is important 
since we would like to know whether the new class of solutions corresponds to 
generic potentials V(<£>) from the particle physics point of view. Therefore, we 
calculate the corresponding inflaton potentials numerically for the scalar case, 
with initial conditions inspired by chaotic inflation. 

This letter is organised as follows. In section II, we briefly remind the reader 
about known analytical solutions. In section III, we present a new family of exact 
solutions and in section IV we study various limits of the parameters of the new 
solutions. 



2 Known analytical solutions 

In the simplest case all coefficients c m vanish and Ust = (Easther's solution 
[BJ). The corresponding functions zs,t are given by ^s,t(^) = Brj + A, where A 
and B are free parameters. For the tensor sector, this gives just the radiation 
dominated Universe. For the scalar sector, the general solution of Eq. (|2]) can 
be written as Hs{v) = G\e n + C%e 7? , where C% and C*2 are two arbitrary 
constants to be determined by the initial conditions (|||). They read C\ — 
and C*2 = — 4y / 7r/pie* fc,?i /v / 2A;. In the large scale limit we have z§ ~ A and the 
spectrum of density perturbations reads: 

e P( = p) 

The spectral index is given by ns = 3, which is in obvious contradiction with 
observations. The analytic form of the corresponding scalar potential V((p) has 
been given in Ref. ||. 

All other cases known so far assume a potential such that L^t = + l)/r/ 2 , 
i.e. c m = if m ^ 2 and c 2 = a(a + 1), where a is a free parameter [ PHI . 
The corresponding functions can be expressed as 2g x(^) = ^4/M a , where A 
is another free parameter. Thus, we have a two-parameter family characterised 
by A and a. For tensor perturbations all these models correspond to power- 
law inflation 0]. For scalar perturbations the subclass a = 1 was studied by 



Starobinsky [fLq| . Another subclass is A = lo\J{a — 1) /a, which corresponds to 
power-law inflation with the scale factor given by a(i]) = l \r]\~ a . The quantity 
Iq has the dimension of a length. In this model inflation occurs if a > 1. In the de 
Sitter case (a = 1), Iq is simply the constant Hubble radius during inflation. Only 
for power-law inflation scalar and tensor perturbations can be solved analytically 
at the same time. For scalar perturbations the general two-parameter family is 
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Figure 1: The generalised Starobinsky potentials for a = 1,2, a p \ = —[3 — 1,5 
(from top to bottom) and A — 1, where a p \ ~ 3.44. The initial conditions are 

7i = (ttpi — l)/«pi ~ 0.7 and if; = ^8nG/3. 

not equivalent to power-law inflation. Of course, a(rj) oc |?7|~ a is a particular 
solution of Zs(rj) = A/\rj\ a , but it is not the general solution. The limit A to zero 
and a close to one reproduces a slow-roll inflation model. Eq. (0) is solved in 
terms of Bessel functions, ^s(v) = (kr})2[CiJ a+ i(kri) + C 2 J_ Q _i (A;^)], where and 
C\ and C*2 are constants fixed from Eq. ([5]). The spectrum may be calculated 
exactly to read 

k 3 Pdk) = ^2 2Q r 2 (a + ~)Ar 2 ^ . (10) 

The corresponding spectral index is ns = 3 — 2a. In particular, we have k 3 P^ = 
lp X /(irA 2 ) and ns = 1 for a = 1. For gravitational waves the same ansatz gives a 
power-law model with spectral index tt-t = 2 — 2a. 

The inflaton potential for the scalar sector is displayed in Fig. (|l]) for various 
values of a. We use initial conditions motivated by the scenario of chaotic infla- 
tion, i.e. the potential and the kinetic energy densities of the inflaton field are 
of the order of the Planck energy density initially. The condition 7; < 1 guaran- 
tees that inflation takes place. For power-law inflation, the potential is given by: 
V(<p) = V^exp[4 1 /7f7(^ - <pi)/m P i]. 
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3 A new analytical solution 



We now turn to the main result of this article. The function zg,T corresponding 
to the effective potential (§) can be written as : 

2Ar 6 ~ 1/2 , / 



where A > and £ = yc 2 + l/4 > 1/2. K 2 ^ is a modified Bessel function of 
order 2£. We now define r as r = 2iki] and k as k = —ici/(2k). c\ defines a 
characteristic scale given by kc = ci/2. Then, the equation of motion (Q) can be 
expressed as: 

dVs, T , r 1 , « , (1/4 -a 



dr 2 



-- + - + 

4 r r 



2 



0S.T = 0. (12) 



The general solution to this equation is given in terms of Whittaker functions 

/x S ,t(t) = CiW^Cr) + C 2 W.^(-r). (13) 

The constants C\ and C2 are fixed by the initial conditions. Using the asymptotic 
behaviour of the Whittaker functions in the small-scale limit, lirn| r |._ >QO W k ^(t) = 
e~ir K , Eq. (9.227) of Ref. pj, we find: 



pikrji+nci / (4k) 

C 1 = T A^l m = , C 2 = 0. (14) 

VZk 



Expressing the Whittaker functions in terms of Kummer functions, Eqs. (9.220) 
of Ref. 1)191 ) we deduce that, if — £ + 1/2 < 0, the growing mode in the large 



scale limit is given by: lim^o W K ^ (r) = T(2£)/[r(l/2 + £ - K)](r)-* +1 / 2 e- T / 2 . 
Using that zs,t(t]) ~ A|?7|2~£ in this limit, it is straightforward to calculate the 
spectrum of density perturbations: 

C ttA 2 2 2 ?" 1 |r(l/2 + £ + ifc c /fc)| 2 ' 1 ; 

This spectrum corresponds to a new exact solution which contains all the 
previously known CclSCS clS particular cases, see the next section. Using Eq. 
(8.328.1) of Ref. 0, we see that for k < k c one has \T(l/2 + £ + ifc c /fc)l 2 « 
27rexp(— 7rkc/k)(kc/k) 2 ^. This implies that in the large-scale limit the spectrum 
is given by: 

7,3 p „ l Pl k 1 ( 2 r 2nkg/k (](]) 

* ^ ~ ttM* (2*0)* ' { } 

if Ci 7^ and cq = 0. We see that it is not analytic in the region of small k. Of 
course, such an infra-red divergence of the spectrum is excluded from observation. 
However, the validity of our description fails for modes that leave the horizon at 
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the Planck epoch, i.e. at the beginning of inflation in the scenario of chaotic 
inflation. This sets a natural infra-red cut-off to the power spectrum (|i"5|). In the 
small-scale limit k ^> kc the spectrum tends to 

A;3p c ~ Jk_2 2 ^- 1 T 2 (Ok 3 - 2i ■ (17) 

7T J\ 

This is the spectrum of the two parameter family studied in the previous section, 
see Eq. fllPl) . The corresponding constant spectral index is given by n$ = 4 — 2£ = 
4 — 2wc 2 + 1/4 and is entirely determined by the coefficient c 2 . 



4 Limiting cases 



4.1 c\ ^ 0, C2 = 0: Coulomb solution 

As already mentioned, the case usually treated in the literature is C\ = 0, c 2 ^ 0. 
On the other hand, the case c x ^ 0, c 2 = has never been studied before. 
Although it is of course just a particular case of the general solution given in the 
previous section, it is worth investigating its properties in some details. For this 
purpose, we restart from the beginning and consider the following function zs t 



z S)T (?7) = 2A y Jc 1 \r l \K 1 \2yJc 1 \r] 

The link with Eq. (|ll]) is obvious since c 2 = corresponds to £ = 1/2. Let us 
now define p and 5 according to p = ir/2, 5 = in. Then the equation of motion 
(H) for Ps,t(v) can b e written as: 

This equation can be solved in terms of Coulomb functions [2(| with I = 0, 
Hs,t(p) = C\Fq{8] p) + C2Go(5; p). The coefficients C\ and C 2 are fixed by the 
initial conditions, see Eq. (||). We find that C\ = ±4i^l Pl e ikm /\/2k, C 2 = 
—iC\. In the large scale limit, the growing mode is given by the irregular Coulomb 
wave function G (S;p) w 1/C (5), where C (S) = ^2n5/ (e 2w5 - 1). Using the 
fact that Ki(x) ~ 1/x when x goes to zero, the spectrum is easily derived. For 
density perturbations, we find: 

fc3 P = ff f e 2«k c /k_ 1 \ (2Q) 



C 2n 2 A 2 k c \ / 
This spectrum is of course nothing but Eq. (|i5|) for £ = 1/2. It is displayed in 
Fig. (Q) for A — 1 and c\ = 0.5. In the regime where k -C kc one recovers Eq. 
(|i"6| ) whereas if k ^> kc this spectrum tends to the Easther's solution k 3 P^ = 
(/p 1 /c 2 )/(7r J 4 2 ) for which = 3 in accordance with Eq. flTT|). The corresponding 
inflaton potentials for various values of C\ are displayed in Fig. (Q). The initial 
conditions are chosen as in Fig. 1 and guarantee that inflation occurs. 
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Figure 2: The scalar power spectrum for the Coulomb case with A = 1 and 
ci = 0.5. 
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Figure 3: The potentials for the Coulomb case for c\ = 0.1, 0.5, 1, 1.5 and A = 1 
with the same initial conditions as before. 
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Figure 4: The scalar power spectrum for the Whittaker case with A 
and Ci = 0.5. 
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4.2 c\ = 0, C2 7^ 0: generalised power-law solution 



Putting c\ = and k = in Eq. (|T5|), one easily checks that the result is 
the generalised power-law spectrum, see Eq. (|10D , with a = £ — 1/2. The 
corresponding spectral index is ns = 4 — 2£. 



4.3 ci ^ 0, c 2 = 2 

Let us finally investigate in more details the case c 2 = 2 which corresponds to 
£ = 3/2. Using Eq. (8.332.1) of Ref. fli~9" |, we find that the exact spectrum for 
density perturbations is given by: 



fc 3p(c 2 =2) 



72 
'PI 



k 



2w 2 A 2 k c 



1 + M 

A; 2 



(21) 



The spectrum is displayed in Fig. (f|) for A = 1 and ci = 0.5. 

In the limit ^c/^ <S 1, the previous spectrum reduces to k 3 P^ C2 ~ 2 ^ = 
(ipi)/ (iiA 2 ), i.e. the spectrum of Starobinsky's solution (ns = 1) as expected. 
The potential of the inflaton is displayed in Fig. (|5|) for various values of c\ with 
the same initial conditions as in Figs. 1 and 3. 

We briefly conclude in recalling that the main result of this article is the 
discovery of a new family of exact solutions for inflationary cosmological per- 
turbations which encompasses all the previously known cases as limiting cases, 
see Table I. This new family could help to shed light on inflationary models 
which do not fulfill the slow-roll conditions. Moreover, it can be used to test 
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Figure 5: The potentials for the Whittaker case for A = 1,C2 = 2 and c± = 
0.1, 0.5, 1, 1.5 with the same initial conditions as before. 



methods of approximation that are valid for models possessing "running" of the 
spectral indices. For this purpose the knowledge of exact solutions allows us to 
avoid a model-by-model analysis, as is necessary for numerical methods. In this 
spirit, the precision of the slow-roll approximation was recently investigated using 
power-law inflation as an exact solution |14| |15] . 
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